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Firmoviscous  and  Anelastic  Properties 
of  Fluids  and  Their  Effects  on  the 
Propagation  of  Compression  Waves 


H.  M.  PAYNTER 


Concept  of  a  Dynamic  Bulk  Modulus 

In  almost  all  dynamical  treatments  of  fluid 
flow,  a  static  relationship  between  the  fluid 
density  and  the  fluid  pressure  is  assumed.  For 
liquids,  this  relationship  is  usually  expressed 
in  terms  of  the  fluid  bulk  modulus,  K,  as  mani¬ 
fested  by: 

=  (1) 


where 

=  fluid  density,  Ib-secV^t^ 

P  =  fluid  pressure,  psf 
K  =  fluid  bulk  modulus,  psf 
On  this  relationship  nearly  all  classical 
developments  for  wave  propagation  in  slightly 
compressible  fluids  have  been  based.  It  is  the 
primary  aim  of  this  paper  to  assume  a  more  gener¬ 
al  dynamical  relationship  between  pressure  and 
density  and  then  cursorily  to  examine  the  result¬ 
ing  effects  on  fluid  phenomena  such  as  acoustic 
wave  propagation. 

From  a  general  macroscopic  and  phenomenolog¬ 
ical  standpoint,  it  is  not  a  priori  reasonable 
to  assume  such  a  purely  static  relationship  be¬ 
tween  pressure  changes  dP,  and  density  changes, 

A;  ,  of  a  fluid  under  dynamic  conditions.  Rather, 
it  is  mo-e  reasonable  initially  to  assume  a  more 
complete  (albeit  linear)  form  such  as 

StfAp  +3iA^  +  a*Aj»  +  -"  =  boAP  +■  b|  AP  +  b2AP+'-- 

(2) 

where  the  dots  indicate  time  derivatives  and  the 
coefficients  are  in  general  undertermlned .  Wils 
relationship  also  may  be  indicated  in  the  equiva¬ 
lent  operational  forms: 


((i.akr'<)Ae  =  Cib»r'')AP 


O) 


with  D  =  d/dt  denoting  the  differential  operator. 

Historically,  only  the  very  simplest  of  the 
infinitude  of  possibilities  for  Equation  (2) 
have  been  considered  analytically  or  have  been 


employed  as  descriptive  models  for  reducing  ex¬ 
perimental  data.  However,  for  a  long  period 
now,  it  has  been  fecognlzed  how  arbitrary  these 
idealizations  have  been;  such  awareness  has  be¬ 
come  particularly  acute  since  the  advent  of  gen¬ 
eralized  continuum  mechanics  and  rheology  (1,  2, 

3).^ 

We  shall  here  consider  only  the  next  approx¬ 
imation  beyond  the  classical;  namely,  that  with 
ao»  ^1»  t>0*  ^1  non-zero,  and  with  all  the  rest 
of  the  a^  and  b;c  assumed  vanishingly  small.  Ihe 
coefficient  aj^  is  directly  related  to  what  has 
been  termed  the  second  coefficient  of  viscosity 

(ii.  1) 

Hiere  is  ample  physical  evidence  to  assume 
that  there  exists  a  time  lag  between  the  appli¬ 
cation  of  pressure  and  the  resulting  change  in 
density.  Using  finite  incremental  changes  this 
phenomenon  might  be  expressed  mathematically  in 
the  assumed  form: 


€ 


r  1  i/h..  1 


r  ocTP»» 
K  L  TP+l 


(4) 


where 

"XT  =  retardation  time,  sec 
D  =  a/dt  =  differential  operation,  sec"^ 
l/K  =  (l/Kj^)  +  (I/K2)  reciprocal  bulk  modulus 
sq  ft  per  lb 

oC  =  a  nondlmenslonal  ratio, 

generally  less  than  unity 
Consequently  a  conjectural  dynamic  or  an¬ 
elastic  bulk  modulus  (which  is  identical  to  the 
static  bulk  modulus  in  the  steady  state)  can 
thus  be  defined: 


To  illustrate  the  physical  significance  of 
the  dynamic  bulk  modulus.  Fig.  1  demonstrates  how 
the  density,  P ,  or  the  pressure,  P,  of  a  fluid 
would  vary  with  time  due  to  the  application  of 
step  change  in  either  quantity. 


Underlined  numbers  in  parentheses  designate 
References  at  end  of  paper. 
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Ca)  S+ep  change  in  density 
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(b)  S+ep  chan^  in  pinessure 

Flg.l  Illustration  of  the  dynamic  bulk 
modulus . 


Effect  on  Wave  Equations; 

Hie  concept  of  a  dynamic  bulk  modulus  can 
now  be  utilized  to  modify  the  evaluation  of  ef¬ 
fective  fluid  capacitance  In  the  derivation  of 
the  wave  equations.  For  the  purpose  of  compari¬ 
son  between  the  pure  wave  propagation  along  a 
pipe  and  the  dispersive  propagation  Involving  a 
dynamic  modulus,  corresponding  variables,  equa¬ 
tions,  and  model  analogies  for  the  two  types  will 
be  carried  in  a  parallel  fashion  as  follows: 

Wondlsoerslve 


dP  -  TiW 
MoaentuBi - ^ 


Elasticity: 
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Fig. 2  Model  analogies  for  a  fluid  element  having 
a  dynamic  bulk  modulus . 


C  =  fluid  capacitance  per  ft  of  pipe,  ft 
t  =  tlm^  sec 

Equations  (6)  and  (7)  can  be  written  In  a 
more  general  form  In  tenns  of  an  Impedance  and 
admittance  operator  that  relate  the  relationships 
between  pressure  and  flow  as  stated  In  Equations 
(8)  and  (9)  for  which  Equations  (6)  and  (7)  are 
merely  special  forms. 


-  -g-  =  z  (o)  ■  w 
_  =  V  C-D)  --p 


Having  written  these  expressions,  the  propa¬ 
gation  operator  and  the  characteristic  Impedance 
Zq  for  a  fluid  column  can  now  be  derived  for  the 
dispersive  case  In  a  similar  manner  to  vdiat  Is 
usually  done  for  the  lossless  case.  Again  for 
the  purpose  of  comparison,  corresponding  relation¬ 
ships  will  be  written  for  the  two  cases: 


Pure  Wave  Frooagatlon 


lapedance 
per  unit 
length 


Admittance 
per  unit 
length 


Z(P)  =  ID 


j  yCv)^ci> 


where.  In  addition  to  the  variables  previously 
defined: 

s  =  distance  along  the  pipe,  ft 

W  =  weight  rate  of  flow,  lb  per  sec 

I  =  fluid  Inertance  per  ft  of  pipe,  sec^/tt^ 


\  r  =  eV^ 

Propagation 

Operator  j  — 'J’p  =  p'*  (*2a) 

Characteristic  I  Zo-Vi/Y 
Impedance  J  — j?*^*^*^ 
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t 
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Pig. 3  Impulse  response  for  an  Infinite 


Dispersive  Propagation 


z(:d)-ip 


(lOb) 


y(:p)=cI^]d 


lb) 
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(12b) 


«rD+i 

where  }.  =  total  column  length,  ft;  T  =  wave  trav¬ 
el  time  along  the  column,  sec.  The  transfer  op¬ 
erator  for  a  pipe  also  may  be  written  In  a  gen¬ 
eralized  matrix  form: 


T. 

W,' 

v. 


cos(t  P  [  P 

^sihli  P  [  eojU  P 


■ 

3 

’ 

wi' 

(U) 


where  subscripts  1  and  2  denote  upstream 
and  downstream  conditions,  respectively. 


Hie  solution  to  Relation  (14)  may  be  obtained  for 
a  given  set  of  conditions  and  fluid  terminations. 
However,  because  of  the  complexity  of  f'  and 
the  general  solution  cannot  be  expressed  In  a 
simple  closed  fonn  as  In  the  case  of  pure  wave 
propagation.  There  are  several  methods  however 
for  approximating  the  solution  for  a  specific 
set  of  conditions. 

Propagation  of  a  Pressure  Impulse  In  a  Seml- 
Inflnlte  Column 

For  example.  If  a  pressure  Impulse  Is  ap¬ 
plied  at  one  end  of  a  seml-lnflnlte  fluid  column. 


Pig. 4  Calculated  response  curves,  (a) 
Impulse  response  I.9O;  (b)  Impulse 
response =  6.33,  15.3;  (c)  step  re¬ 
sponses  . 


the  shape  of  this  Impulse  at  any  point  can  be 
approximated  by  certain  distribution  functions 
as  It  travels  along  the  column. 

Since  If  such  a  seml-lnflnlte  column  the 
pressure  at  any  point  is  given  by: 

FCs.e)  =  e'’ •  (15) 

This  result  follows  from  Equations  (12)  with 

HoW(Sy<)  =  (16) 

which  Is  true  for  the  case  In  question.  Follow¬ 
ing  the  methods  of  one  of  the  authors  J_) , 
the  operator  c  may  be  expanded  as  follows: 


3 


In  Equation  (17)  the  time  constants  have  the 
subsequent  significance 

Tj,  =  T  =  mean  delay,  sec 

Tg  =  J^/^=dlsperslon  time,  sec 

Hanln  (2)  and  Vaughn  (8)  have  demonstrated 
for  the  case  '  =  0  that  an  Impulse  at  the  ori¬ 

gin  gets  Increasingly  dispersed  and  attenuated 
as  It  propagates.  Following  suggestions  of  the 
present  authors,  Vaughn  showed  that  Its  shape 
can  be  predicted  accurately  by  the  delayed  chi- 
square  distribution  function  (II)  and  that  the 
limiting  form  for  very  early  development  can  be 
predicted  by  the  complementary  error  function 
(I).  For  large  Intervals  of  time  the  pulse  ap¬ 
proaches  the  normal  distribution  function  (III). 

Figs.  5  and  4  Illustrate  these  shapes 
graphically,  while  the  Appendix  presents  the  cor¬ 
responding  mathematical  formulas,  for  the  impiilse 
and  step  response  of  a  seml-lnflnite  fluid  medi¬ 
um. 

However,  the  two  derived  constants,  T  and 
Tg,  adequately  define  behavior  for  all  types  of 
disturbances  and  boundary  conditions.  This  fact 
permits  determination  of  the  fluid  properties 
from  any  adequate  and  consistent  observational 
data. 


Conclusions 

All  experimental  data  known  to  the  authors 
reinforce  the  belief  that  such  a  compressive 
model  as  that  proposed  In  the  foregoing  Is  of 
the  minimum  complexity  which  will  still  explain 
the  scattering  action  readily  apparent  In  ob¬ 
served  behavior  of  fluids.  Thus  we  strongly 
urge  that  a  concerted  effort  be  made  to  reduce 
such  observations  to  the  point  where  the  anelas- 
tlc  properties  of  the  more  common  fluids  can  be 
established  firmly. 
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APPENDIX 

Derivation  of  Dynamic  Response  Equations 

Purely  Flrmovlscous  Propagation  ( ''■  =  0). 
For  this  case  it  Is  possible  to  simplify  analy¬ 
sis  by  first  normalizing,  or  nondlmenslonallz- 
Ing,  all  variables  and  parameters  through  Intro¬ 
ducing  Into  Equation  (15)  the  substitutions: 

NORMALIZED  DISTANCE:  5C  S  (T/^)  •  X 
NORMALIZED  TIME:  ^  =  ( ’/^  )  •  ^  (A) 

NORMALIZED  DERIVATIVE:  <£>  =  C  ® 


Thus  we  may  now  treat  the  generalized  operator; 


e''  =F(^,£>)  =  6 


■‘f , 

The  condition  03  corresponds  to 

the  instantaneous  response,  t  -►  0,  and  yields: 

=Foo  (0C)  =  e-^^  (J 
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nils  operator  applied  to  a  unit  Impulse  or 
unit  step  results  In  the  consequent  Impulse  re¬ 
sponse  and  step  response.  In  the  form: 

IMPULSE  RESPONSE;  fSC,  ^  '^7' 
STEP  RESPONSE:  Ko  (X,V~)  =  1  “  er/  (^/yfz^ 
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Tables  of  the  functions  Indicated  are  readi¬ 
ly  available  In  standard  handbooks,  such  as  that 
of  Burington  and  May  (2^) . 

For  the  opposite  conditions,  when  the  dis¬ 
turbance  has  traveled  a  great  distance  from  the 
origin,  then  t-»  co,  iD— >  0,  and  we  may  consider 
the  limiting  operator; 


Between  these  two  extremes  many  models  would 
be  possible.  That  used  by  Vaughan  (2.),  as  men¬ 
tioned  earlier,  is  a  “delayed  chi-square  distri¬ 
bution,"  which  is  an  approximation  in  the  form: 

(F) 


corresponding  to  the  Impulse  and  step  responses: 

IMPULSE  RESPONSE:  (OC,  Vj  = 


VjtTt^ 
2- 


(E) 


STEP  RESPONSE: 


with  Tj,  T2,  normalized  time  constants  and  m  the 
number  of  degrees  of  freedom  of  the  chi-square 
distribution.  An  excellent  fit  is  obtained  by 
matching  the  first  three  Impulse  moments,  which 
yields  the  conditions: 

t;  =  =-  Q.xwoc  (G) 

T1  =  ^  =  1-125 

w  =  =  1.580  2)C 


Thus  these  results  correspond  to  the  well-known 
Gaussian  distribution  or  normal  probability  dis¬ 
tribution,  with  mean  =  variance  =x-  This  func¬ 
tion  is  tabulated  in  many  handbooks,  such  as 
that  previously  cited  (^) . 


The  corresponding  values  of  Impulse  response, 
^l/m*  step  response  ure  readily  found 

from  any  of  the  available  tables  of  the  chi- 
square  distribution  such  as  those  previously 
cited  (IB) . 
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